Pair Structure and the Pairing Interaction in a Bilayer Hubbard model 
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The bilayer Hubbard model with an intra-layer hopping t and an inter-layer hopping t± provides 
an interesting testing ground for several aspects of what has been called unconventional super- 
conductivity. One can study the type of pair structures which arise when there are multiple Fermi 
surfaces. One can also examine the pairing for a system in which the structure of the spin-fluctuation 
spectral weight can be changed. Using a dynamic cluster quantum Monte Carlo approximation, we 
find that near half-filling, if the splitting between the bonding and anti-bonding bands t±/t is small, 
the gap has Bi g (d x 2_ y 2-wave) symmetry but when the splitting becomes larger, A\ g (s -wave) 
pairing is favored. We also find that in the a pairing region, the pairing is driven by inter-layer 
spin fluctuations and that T c is enhanced. 
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The bilayer Hubbard model provides a model system 
for which one can study the relationship between the 
structure of the Fermi surface, the spin-fluctuation spec- 
trum and superconductivity. By varying the relative 
strength of the inter-layer one electron hopping t± to 
the near-neighbor intra-layer hopping i, one can alter the 
size and shape of the bonding and antibonding Fermi sur- 
faces, change the momentum and frequency structure of 
the spin-fluctuation spectral weight, and move away from 
the Mott region. Previous Monte Carlo calculations [lj- 
6] for a doped bilayer found evidence for an attractive 
pairing interaction in both the d x 2_ y 2 (cosk x — cosfcj,) 
and cos k z channels. The latter was called a d z chan- 
nel and corresponded to a pairfield function which had 
different signs on the bonding and antibonding Fermi sur- 
faces. This is an Ai g gap and here we will refer to it as 
an s ± gap. The sign problem associated with the pre- 
vious determinantal quantum Monte Carlo calculations 
prevented one from numerically exploring the low tem- 
perature properties of the doped bilayer. More recently, 
a functional renormalization group study [7] found that 
as t±/t increased, s pairing was favored over d x 2_ v 2 
pairing. Here using a dynamic cluster approximation^] 
(DCA) we will examine what the bilayer Hubbard model 
tells us about the pairing mechanism and the search for 
higher T c materials. 

The Hamiltonian for the bilayer Hubbard model shown 
in Fig. [T]can be written as 

H = -t ^ ( c ]m/™ + h.C.) - t±^2{cf lo .Ci 2 <r + h.C.) 
{ij)ma i<J 

- ft/ J n im <7 + U 22 n irn-fniml (1) 
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Here t and t± are intra- and inter-layer near neighbor 
hopping parameters, /i is the chemical potential and U 
is an onsite Coulomb interaction. The indices i and j 
run over the sites in both the m = 2 (upper) and m = 1 
(lower) layers and (ij) implies that only a single (ij) near 
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FIG. 1: a) The bilayer Hubbard lattice with near neighbor 
intra-layer t and inter-layer t± hopping parameters and an 
onsite U Coulomb interaction. The bonding (fc z = 0) and 
antibonding (k z — it) Fermi surfaces for t± — 0.5 (b) and 
2.0 (c) at a filling (n) — 0.95. A d x 2_ y 2 gap structure is 
illustrated for the t±/t = 0.5 Fermi surface and an s gap 
structure is shown for t±/t = 2.0. 



neighbor hopping is included in the sum. In the following 
we will measure energies in terms of t and set U = 6. The 
bonding (fc 2 = 0) and anti-bonding (k z — n) bands are 
given by 



e(k) = — 2i(cos k x + cos k y ) ± t± cos k z 



(2) 



The k z = bonding and k z — ir anti-bonding Fermi 
surfaces (FS) of the non-interacting system are shown in 
Fig. [1] for a filling (n) = 0.95 and two values of t±/t. 
There is an anti-bonding (k z = n) electron-like FS that 
forms around the (k x = 0, k y = 0) T point of the 2D Bril- 
louin zone and a bonding (k z — 0) hole-like FS around 
the (it, 7t) point. For (n) = 0.95, the topology of the 
non-interacting FS changes from two electron-like FS's 
to one electron and one hole- like FS when t±/t exceeds 
a critical value of order 0.1. Here we will focus on the 
behavior of the bilayer system for t±/t greater than this 
value, although in some plots we will give the t± = re- 
sult. As t±/t increases further, the FS's shown in Fig. 1 



shrink and for (n) = 1 the non-interacting system be- 
comes a band insulator when t± > 4i. Weak coupling 
fluctuation-exchange (FLEX) [l| as well as phenomeno- 
logical 9] spin-fluctuation calculations find superconduc- 
tivity for the doped system with a d x 2_ y 2-\ike gap for 
t±/t = 0.5 and an s ± -like gap for t±/t = 2.0, as schemat- 
ically illustrated in Fig. [1] This two- Fermi-surface model 
can be seen as an analog of the multi-Fermi-surface Fe- 
pnictide models. 

In the following, we will use the DCA to study the 
Hubbard bilayer with U = 6 and (n) = 0.95. The DCA 
[g, [l0( maps the bulk lattice problem onto an effective 
periodic cluster embedded in a self-consistent dynamic 
mean-field that is designed to represent the remaining 
degrees of freedom. The DCA calculations were carried 
out on an N = 2 x (L x L) cluster with L — 4 and the 
effective cluster problem was solved using a Hirsch-Fye 
quantum Monte Carlo algorithm [10j. 

The magnetic susceptibility is given by 



X (q)= / dr<m-(r)m+(0)) 



Here 



and 



m q = 7^E c k+ q t c H 
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Fig. [2] shows plots of x(tt, tt, 0) and \ 1 (0) versus t± at a 
temperature T = 0.4. For the undoped system at T = 0, 
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FIG. 2: The spin susceptibility x( 7r i 7r )0) and X _1 (0) for 
(n) = 0.95, (7 = 6 and T = 0.4 are plotted versus t±. As 
£_L exceeds 2, the planar AF response is suppressed and sin- 
glet correlations between sites on opposite sides of the bilayer 
become dominant. 

the AF order vanishes for t± > 2.0 and the ground state 
is a disordered valence bond (VB) phase with a spin 



gap [6]. Here, for the doped system, one sees that for 
t± ~ 2 the near neighbor in plane (tt, it, 0) response is 
suppressed and as t± increases further, the low energy 
interlayer spin-fluctuations become gapped as the inter- 
layer valence-bond singlets form. This crossover is also 
clearly seen in the behavior of the inverse spin suscepti- 
bility x ~Hq = 0). 

As noted, weak coupling calculations have found both 
d x 2_ y 2 (Big) and s^ 1 (A lg ) pairing correlations for the 
Hubbard bilayer model. Here, in order to probe the su- 
perconducting response we have used the DCA to calcu- 
late the pairfield susceptibilities. 



P«(T) 



dr<A Q (r)AL») 
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associated with each of these symmetries. For the d x 2_ y 2- 
wave we have taken 



A l 2 -y 2 = -^^2(cos k x - cos k y )cl t c 
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and for the s^ 1 case [11 1 



A 1 



1 



vr 



cosfc z cL.c 



2 °fet -k]. 



(7) 



(8) 



Here, k = (k Xl k yi k z ) with k z — and 7r. In Fig. [3l 
results for these pairfield susceptibilities are shown for 
different values of t± versus temperature. For t±_ = 0.5, 
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FIG. 3: a) The s and d-wave pairfield susceptibilities P a 
versus temperature for different values of the inter-layer hop- 
ping t±. As t± increases, the leading pairfield susceptibility 
changes from d-wave to s and at larger values of t± the su- 
perconducting pairfield susceptibility is suppressed as inter- 
layer valence bonds form. 

one sees in Fig. [3^ that the d x 2_ y 2-waMe susceptibility 
is rising the most rapidly at low temperatures followed 
by the s -wave response. In Fig. [2b, for t± = 1.0 the 



d-wave response is slightly larger than the s 1 * 1 case but 
they are closer, while for t± = 2.0 the s^ is the leading 
pair susceptibility. This result is similar to that found 
in Ref. [7]. The final panel, Fig. [3Ji, shows the pairficld 
susceptibilities for t± = 3. For this value of t±, the T = 
undoped system is well into the spin-gapped VB phase 
and one sees that the pairficld response of the doped 
system is significantly weakened. 
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4: (a) The strength of the pairing interaction V a for a = 
d and s calculated from Eq. ((9]) versus t±. The integrated 
spectral weights I v , Eq. (|10[) . for the intra- (y — x) and inter- 
[y = z) layer near-neighbor spin- fluctuations versus t±. 



The pairing interaction is given by the irreducible 
particle-particle scattering vertex T(K, K') which de- 
scribes the scattering of a singlet pair from the state 
(K,-K) to (K',-K') with K = (K,iw„). Given this 
vertex and the dressed single-particle Green's functions, 
one can solve the Bethe-Salpeter equation to find the 
pairing eigenvalues Aq, and eigenfunctions tp a (K) [l2j. 
For the bilayer Hubbard model, the leading pairing eigen- 
value at small values of t± corresponds to the B\ g 
(d x 2_ y i) state and, as t± exceeds 1, switches to the Ai g 
(s ) state. A measure of the pairing strength V a for a 
given pairing state is 



V a 



l / N Ek,k- <A*( K > ^ r ) r ( K > *T, K ' *T)<p a {K', ttT) 

(9) 
A plot of V a versus t± is shown on the left hand side of 
Fig. 4. On the right hand side of this figure we have plot- 
ted the integrated spin-fluctuation spectral weights for 
the intra- and inter-layer near-neighbor spin fluctuations 
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Here K„ = K^ for the intra-layer fluctuations and K 2 
for the inter-layer fluctuations, and x c (K,cj) is the spin- 
susceptibility calculated on the 2 x (4 x 4) cluster. As 
t± increases, the dominant spin-fluctuations change from 



intra-layer to inter-layer. The inter-layer spin fluctu- 
ations give rise to the scattering of pairs between the 
bonding and anti-bonding Fermi surfaces and lead to s 
pairing. We believe that the similarity of the t± depen- 
dence of I v and the pairing strengths V a provides evi- 
dence linking the structure of the spin fluctuation spec- 
tral weight to the pairing mechanism and the resulting 
gap structure in this model. This spin-fluctuation pic- 
ture differs from an alternate strong coupling scenario in 
which U is assumed larger than the bandwidth and su- 
perconductivity is viewed as arisin g fr om the melting of 
a valence bond phase with doping [13J]. We believe that 
the spin-fluctuation picture provides a more appropriate 
framework for U of order the bandwidth, which corre- 
sponds to the parameter region in which the pairing is 
strongest. 

Although we cannot carry out the finite size scaling 
analysis one would need to determine T c , we have esti- 
mated T c values from the temperature at which the DCA 
results for P a (T) diverge. We find that the maximum T c 
for the s^ phase is significantly larger than that for the 
d-wave phase even though, as seen in Fig. 4, the maxi- 
mum value of Vd exceeds that of V s ±. This is because, 
at (n) — 0.95, the t±_ —¥ system is near the Mott state 
and the quasi-particle weight is strongly suppressed. At 
larger values of t±, where V s ± peaks, the bilayer sys- 
tem with an intermediate U = 6 has moved into a semi- 
metallic phase which has only a moderate renormaliza- 
tion of the quasi-particles. It is interesting to note that 
the maximum T c for the s phase with (n) = 0.95 is more 
than a factor of two larger than the T c obtained from a 
corresponding DCA estimate of the single-layer case at 
optimal doping. In this case, the T c of the single-layer 
system has been increased by doping away from half- 
filling so as to suppress the quasi-particle Mott renor- 
malization. However, this reduces the strength of the 
pairing interaction. In the bilayer case, one can move 
away from the Mott region by increasing £j_ and at the 
same time still achieve a large interaction strength as 
shown in Fig. 4. Thus we conclude that multi Fermi sur- 
face materials can offer a possible pathway to higher T c 
superconductivity. 

To conclude, in this paper we have reported DCA re- 
sults for the pair structure and pairing interaction for a 
bilayer Hubbard model. By increasing the relative size of 
the inter-layer hopping t± to the intra-layer hopping t at 
a fixed doping, we have shown that one can tune the sys- 
tem from a d x 2_ y 2 (Bi ff ) superconductor to an s ± (A lg ) 
superconductor. We then examined the strength of the 
pairing interactions in these channels as t±/t was varied 
and showed that they were correlated with the relative 
strengths of the near-neighbor intra-layer to inter-layer 
spin fluctuations respectively. When the inter-layer spin 
fluctuations become dominant, an s superconducting 
phase is favored. The ability to change the nature of the 
superconducting state by varying a one-electron param- 



eter in the Hamiltonian represents a useful approach for 
studying the relationship between the pair structure and 
the underlying pairing interaction as well as identifying 
the basic correlations which are responsible for pairing. 
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